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ABSTRACT 

Combining a theoretical model of mass accretion onto a galactic center with a high-resolution N- 
body/SPH simulation, we investigate the formation of an intermediate massive black hole (IMBH) during 
the hierarchical formation of a small spiral galaxy (with a total mass of 1O 1O M0) in the high-z universe. 
We found that the rate of average mass accretion to the nucleus due to the radiation drag exerted by 
newly formed stars in the forming galaxy is w 10~ 5 M Q yr~ 1 . As a result of this accretion, an IMBH 
with « 10 4 M Q can be formed in the center of the spiral galaxy at z ~ 4. We found that a central 
BH coevolves with the dark matter halo from z^l5toz^2. The mass ratio of the BH to the dark 
matter halo is nearly constant « (1 — 3) X 10 -6 from z ~ 10 to z ~ 2. This is because that change 
in the dark matter potential enhances star formation in the central part of the galaxy, and as a result 
the BH evolves due to mass accretion via the radiation drag. Therefore, our model naturally predicts a 
correlation between massive BHs and dark matter halos. Moreover, it is found that the hnal BH-to-bulge 
mass ratio («5x 1CU 5 ) in a small spiral galaxy at high-z is much smaller than that in the large galaxies 
(« 10 -3 ). Our results also suggest that the scatter in the observed scaling relations between the bulge 
mass and black hole mass are caused by a time lag between BH growth and growth of bulge. We also 
predict that the X-ray luminosity of AGN is positively correlated with the CO luminosity in the central 
region. By comparing our results with the properties of Lyman break galaxies (LBGs), it is predicted 
that some LBGs have massive BHs of sa 10 6 — 1O 7 M . 

Subject headings: black hole physics - galaxies: nuclei, starburst - hydrodynamics - radiation 
mechanisms: general - method: numerical 



1. INTRODUCTION 

Recent compilation of the kinematical data on galactic 
centers has revealed that a central "massive dark object" 
(MDO), which is a candidate for a supermassive black hole 
(BH), tightly correlates with the mass of a galactic bulge; 
the BH-to-bulge mass ratio is w 0.001 as a median value 
(e.g., Kormendy & Richstone 1995; Magorrian et al. 1998; 
Merritt & Ferrarese 2001; McLure & Dunlop 2002; Mar- 
coni & Hunt 2003). There have been a number of theoret- 
ical efforts to clarify the origin of this relation (e.g., Silk 
& Rees 1998: Ostriker 2000; Adams, Graff & Richstone 
2001). However, little has been elucidated regarding the 
physics on the angular momentum transfer in a spheroidal 
system (a bulge), which is inevitable for formation of BHs. 
Recently, Ferrarese (2002) and Baes et al. (2003) have ar- 
gued that the BH mass in spiral galaxies is related to the 
dark matter halo mass. This correlation suggests that for- 
mation of supermassive BHs is physically connected not 
only with formation of galactic bulges, but also with as- 
sembly processes of dark matter halos in galaxy formation. 
Since merging of proto-galaxies triggers active star forma- 



tion, a physical link between star formation and mass ac- 
cretion toward the central BH is expected. 

Umemura (2001) has considered the effects of radiation 
drag as a mechanism for removing the angular momen- 
tum of the gas in the active galactic nuclei. The radia- 
tion drag in the solar system is known as the Poynting- 
Robertson effect. Note that, in the early universe, Comp- 
ton drag force has a similar effect on formation of massive 
BHs (Umemura, Loeb, & Turner 1997). The rate of an- 
gular momentum loss due to radiation drag is given by 
dhxJ/dt ~ —XdE/c, where J is the total angular momen- 
tum of the gaseous component, E is the energy density of 
the uniform spheroidal system, and Xd is the mass extinc- 
tion coefficient which is given by Xd — ^d^d/Pgas with the 
number density na, cross-section Cd and gas density p gas - 
The exact expressions for the radiation drag are given in 
the literature (e.g., Umemura, Fukue, & Mineshige 1997; 
Fukue, Umemura, & Mineshige 1997). 

In an optically thin regime, dhxJ/dt ~ — (rL st /c 2 M gas ), 
where r is the total optical depth of the system, L* is 
the total luminosity of the spheroidal system, and Af gas is 
the total mass of gas. In an optically thick regime, the 
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radiation drag efficiency is saturated due to conservation 
of the photon number (Tsuribe & Umemura f997). Thus, 
an expression of the angular momentum loss rate suitable 
for both regimes can be d\n J/dt ~ — (L*/c 2 M g )(l — e~ T ). 
The mass accretion rate is therefore M — —M g d\n J/dt = 
(X»/c 2 )(l — e~ r ). In an optically thick regime, this gives 
simply M = L*/c 2 (Umemura 2001). Thus, the to- 
tal mass accreted onto the MDO, Mmdo, is maximally 
Mmdo — / L*/c 2 dt. For more realistic cases, we should 
take into account inhomogeneity of the interstellar matter 
(ISM). In active star- forming galaxies, the ISM is observed 
to be highly inhomogencous (Sanders et al. 1988; Gor- 
don, Calzetti & Witt 1997). In addition, high resolution 
three-dimensional hydrodynamic simulations have shown 
that multiple supernovae (SNe) in a galactic center form a 
quasi-stable inhomogeneous torus around a supcrmassive 
black hole (Wada & Norman 2002). In such inhomoge- 
neous ISM, optically thin surface layers of optically thick 
clumpy clouds lose their angular momentum due to radia- 
tion drag, and eventually they accrete toward the galactic 
center (Sato et al. 2004). Kawakatu & Umemura (2002) 
have shown that the inhomogeneity of ISM plays an impor- 
tant role in the radiation drag attaining maximal efficiency. 
Based on the radiation drag model, Kawakatu, Umemura, 
& Mori (2003) predict that a mass ratio between the black 
hole mass and the bulge mass, MsH/^buige — 0.001, which 
is determined by the energy conversion efficiency of nu- 
clear fusion from hydrogen to helium, i.e., 0.007 (Umemura 
2001). In these previous studies, galaxies are treated as a 
one-zone model, therefore the growth of supermassive BHs 
has not been revealed in a more realistic situation, namely 
that of the hierarchical formation of galaxies. Granato et 
al. (2004) presented a semi-analytic modeling of the early 
evolution of massive spheroidal galaxies and AGNs within 
the dark matter halo, including the angular momentum 
transfer via radiation drag. They claimed the feedback 
from supernovae and from AGNs determines the relation 
between the BH mass, the bulge mass and the dark matter 
halo mass (see also Bukert & Silk 2001). 

Di Matteo et al. (2003) followed the evolution of the 
gas, stars, and the dark matter in forming galaxies, and 
they found that the observed BH mass-to-stellar velocity 
dispersion of a bulge is reproduced if the gas mass in the 
bulge is linearly proportional to the black hole mass. How- 
ever, it is impossible to resolve the structure of the central 
sub-kpc region or the bulge component of host galaxies be- 
cause of the limitation on their numerical resolution (mass 
resolutions are 10 7 — 10 8 M Q and gravitational softening 
lengths is 4-9 kpc). Recently, using high-resolution cos- 
mological A-body/SPH simulations with 2 x 10 6 particles 
(one SPH particle has 10 3 M Q and a gravitational soften- 
ing length of ~ 50 pc), Saitoh & Wada (2004) investigated 
stellar and gaseous cores on a sub-kpc scale during the 
hierarchical formation of a small spiral galaxy (with a to- 
tal mass of 10 10 M Q ), and found that the galactic core (< 
several 100 pc) coevolves with the galactic dark matter 
halo of ~10 kpc scale. The rapid increase of the gas mass 
and star formation rate (SFR) in the central sub-kpc re- 
gion may cause further mass accretion to the nucleus, due 
to, for example, a turbulent viscosity (Wada, Meurer, & 
Norman 2002), a gas drag and dynamical friction in dense 
stellar clusters (Norman & Scoville 1988), the radiation 



drag originating in the nuclear starburst (e.g., Umemura 
2001), or the BHs-BHs merger (e.g., Haehnelt 2004). 

Here we focus on the radiation drag as one of the pos- 
sible processes of mass accretion onto a BH during hi- 
erarchical formation of a galaxy. The observed AGN- 
starburst connection in nearby galaxies (Heckman et al. 
1989; Kauffman et al. 2003; Imanishi & Wada 2004; 
Jahnke et al. 2004) suggests that star formation plays 
an important role in the mass accretion. We expect a cor- 
relation between BHs and bulges as a natural result of 
galaxy formation, if the radiation drag works. In this pa- 
per, we quantitatively estimate evolution of a black hole 
mass in a forming galaxy, combining A-body/SPH simu- 
lations of cosmological galaxy formation done by Saitoh & 
Wada (2004) with an analytic model of angular momen- 
tum transfer due to the radiation drag. 

This paper is organized as follows: In Section 2, we 
briefly describe the simulation of galaxy formation. Our 
model for the growth of a massive BH via the radiation- 
hydrodynamic process is also explained. Based on this 
model, in Section 3, we show the history of accretion to 
a BH in a spiral galaxy. We also discuss the mutual rela- 
tionships between a massive BH, a galactic bulge, and a 
dark matter halo. Finally, we discuss correlation between 
AGN activities and the properties of bulges. In Section 
4, we compare our predictions with observational scaling 
relations, and we discuss the Lyman break galaxies, as can- 
didates for the small spiral galaxies that we demonstrate 
here. Section 5 is devoted to conclusions. 

2. MODELS 

2.1. Simulations of galaxy formation 

The numerical simulations used here are based on Saitoh 
& Wada (2004) and Saitoh et al. (in preparation). We 
model the formation and evolution of galaxies in the CDM 
universe, adopting a top-hat initial condition with an open 
boundary for a single galactic halo (Mhaio ~ 10 10 M Q ). 
The cosmological parameters in our model are fi = 1.0, 
n x = 0.0, n h = 0.1, h = tfo/km/s/Mpc = 0.5, and 
as = 0.63. The collapse epoch of the halo is set at z c <~ 3 
and its spin parameter is 0.05 (Barns & Efstathiou 1987; 
Heavens & Peacock 1988). 

Since the total mass of the object in our simulation is 
small (10 10 M Q ), evolution of the system and therefore the 
conclusion in this paper do no depend on the employed 
cosmology. The collapse epoch of the object in our sim- 
ulation is z c ~ 3, for which the evolution is not strongly 
affected by the A term. This is in contrast to much larger 
systems, such as clusters of galaxies. 

The number of baryon (SPH) and dark matter parti- 
cles in the spherical region is Asph = Adm = 1005600. 
The mass resolutions of baryon (gas and stars converted 
from the gas) and DM particles are 1.1 x 10 3 M Q and 
1.0 x 10 4 M Q , respectively. The gravitational softening 
lengths are 52 pc for baryon particles and 108 pc for DM 
particles. The initial distribution of the particles is gen- 
erated by COSMICS (Bertschingcr 2001). We discuss the 
evolution of galaxies until z = 2, because the assembly 
history for the boundary conditions would not be realistic 
much later than the collapse epoch, z c . But, z ~ z c the 
assembly of the simulated galaxy finishes, so that the mass 
of the galaxy at z — would be equal to that at z = 2. 
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bulge scale. The N[ nt is defined by iVi nt = n c 7rr 2 rbui g c = 
(3/4)iV c (r c /rbuigo) 2 , where n c = N c /§7rrg ulge is the num- 
ber density of gas clouds. In this paper, we assume that 
the cloud covering factor is order unity, i.e. N- mt w O(l), 
according to the previous analysis. A different level of ISM 
dumpiness can reduce the radiation drag efficiency by a 
factor of 2 (Kawakatu & Umemura 2002). We also confirm 
that the optical depth in a clumpy media is an order of 
unity, using a 3-D hydrodynamic simulations (see Fig. 1 
in Wada & Norman 2002). Even if the system is extremely 
gas-rich (M gas = 10 8 M Q in the central 100 pc), we found 
that the optical depth for the disk plane from various di- 
rections is distributed between 0.5 and 1.2, assuming the 
same gas/dust ratio and dust opacity in the analysis here. 



Finally, using M { 



gas 



N c m gas the total optical depth of 



the bulge can be re-written to be 

TbulgeW = f(t)N int ~ 



3xd M gas (t) 



An r? 



The numerical technique we employ to represent 
the evolution of galaxies is a standard hybrid N- 
body/hydrodynamic code for galaxy formation. The code 
includes both the radiative cooling and star formation. 
However, the dynamical and radiative feedback processes 
from star formation and supernova explosions are not ex- 
plicitly taken into account. The length of the interaction 
list of each SPH particle is N^b = 50. In the SPH sim- 
ulations, the Jeans instability can be resolved correctly 
for masses larger than 27Vnbwsph (Bate & Burkert 1997), 
where tosph is mass of an SPH particle. In the simulation, 
we can resolve gravitational instability of a cloud whose 
mass is larger than 1.1 x 10 5 M Q . In order to model the 
multiphase nature of the interstellar medium (e.g., Wada 
& Norman 2001), we solve the energy equation with the 
radiative cooling under 10 4 K and the inverse Compton 
cooling. We assume that the gas has a primordial abun- 
dance of X=0.76 and Y=0.24 and we assume an ideal gas 
with 7 = 5/3. The mean molecular weight of gas, /i, is set 
to 0.59. 

The star formation algorithm is similar to the one by 
Katz (1992). If an SPH particle satisfies all the fol- 
lowing conditions: (1) the regions are in virialized halos 
(psph > 200/9bg)j where pec is the background density, 
(2) low temperature (T< 3 x 10 4 K), and (3) collapsing 
regions (V • v < 0), then it is converted into a collision- 
less star particle inheriting the velocity and the mass of 
the gas particle. The local star formation rate, SFR, is 
assumed to be SFR=c*mspH/Tff with = 1/30, where 
Tfi = 1/y/GpspH is a free-fall time. 

2.2. Radiation drag model 

Next, we model the BH growth, based on the radia- 
tion drag-driven mass accretion. Here, we suppose a two- 
component system: an inhomogeneous ISM with disk-like 
geometry embedded in a spheroidal stellar bulge. In fact, 
the ISM in active star-forming galaxies is highly inhomo- 
geneous (Sanders et al. 1988; Gordon, Calzetti & Witt 
1997). Wada & Norman (2002) also showed that the mul- 
tiple SNc in a galactic center from a quasistable inhomo- 
geneous torus around a supermassivc BH. 

We assume that stars and N c (= 10 4 ) identical clouds 
are randomly distributed in a system of the bulge. 2 The 
optical depth of a gas cloud is f = Xd/Ogas^c — Xdm gas /r 2 , 
where p gas , w gas , and r gas are the densitiy, mass and the 
size of a cloud. The mass extinction coefficient Xd is given 
by Xd = «dCd/pgas with the number density rid, cross- 
section <7d and gas density p gas . In this paper, we assume 
Xd = 300cm 2 g- 1 (ad/O.l/imr 1 (p s /gc m - 3 ) (Z/O.3Z ), 
where ad is the grain radius, p s is the density of solid 
material density within the grain (e.g., Spitzer 1978), and 
Z is the metallicity of gas, which are fixed at ad = O.l/i m, 
p s = lgcm~ 3 and Z=0.3Z Q . 3 The total optical_depth of 
the bulge T bu igo is given by T bu i g0 = fN int , where N int is the 
average number of clouds intersected by a light ray over a 

2 It is noted that simulations with a three times larger number of clouds did not lead to any fundamental difference for final BH mass, although 
at least 10 4 clouds are necessary to treat the radiation transfer effect properly in clumpy ISM. The total optical depth, Tb u lge> is not also 
significantly affected by changing the cloud size, r c . 

3 We should keep in mind that recent observations suggest that the metallicity of the gas in the AGNs can be super-solar Z > Zq (e.g., Ohta 
et al. 1996; Dietrich & Wilhclm-Erkcns 2000; Maiolino ct al. 2003). If this is the case, the optical depth of a gas cloud can be enhanced by a 
factor of 3-4. 

4 In the present paper, we identify a 'bulge' as a spheroidal star-forming region where the average number density of the gas, «h, is larger than 
0.1 cm -3 in a spiral galaxy. This criterion corresponds to the density criterion of the star-forming region. 



(tV (1) 

' bulge W 

where rb u \ ge (t) and M gas (t) are the size and the gas mass 
of the bulge. 4 

The radiation drag, which drives the mass accretion, 
originates in the relativistic effect in absorption and subse- 
quent re-emission of the radiation. This effect is naturally 
involved in relativistic radiation hydrodynamic equations 
(Umemura, Fukue, & Mineshige 1997; Fukuc, Umemura, 
& Mineshige 1997). The angular momentum transfer in 
radiation hydrodynamics is given by the azimuthal equa- 
tion of motion in cylindrical coordinates, 

l_d{rv^ ^Xd _ {E + P H> )v l (2 ) 

r at c 

where E is the radiation energy density, is the radiation 
flux, is the radiation stress tensor. By solving radia- 
tive transfer with including dust opacity, we evaluate the 
radiative quantities, E, F^, and P^, and thereby obtain 
the total angular momentum loss rate. Then, we can esti- 
mate the mass accretion rate of the dusty ISM accreted on 
to a central massive dark object, Mmdo, by using the rela- 



tion, M drag /M g 



-J I J, where J and M gas are the total 



angular momentum and gas of ISM. In the optically-thick 
regime of the radiation drag, —M gas J/J = L\j U \ gc (t) / c 2 , 
where Lbuige(i) is the total luminosity of the bulge. The ra- 
diation drag efficiency depends on the optical depth Tb u i g e 
in proportion to (1 — e ^ Tbuls °^) (Umemura 2001). Thus, 
the mass accretion rate via the radiation drag is 

^^^^^(l-e-W*)), (3 ) 

where ibuigo(^) and Tb u i ge (i) are the total luminosity 
and the time-dependent total optical depth of the bulge. 
Kawakatu & Umemura (2002) found that the efficiency 
?7drag is maximally 0.34 in the optically thick regime. 

The radiation energy emitted by a main sequence star is 
0.14e to the rest mass energy of the star, where e is the en- 
ergy conversion efficiency of the nuclear fusion from hydro- 
gen to helium, which is 0.007. Thus, the luminosity of the 
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bulge at optical and UV bands is simply approximated by 
£buige(*) - 0.14eM b uigo(t)c 2 , where M bu igc(£) is the SFR 
in the bulge. Here, we employ a stellar initial mass func- 
tion (IMF) such as <j> — dn/d\ogm* = A(TO*/M ) _a for 
a mass range of [mi, m u ], where to*, mi, and m u are the 
stellar mass, the lower mass, and the upper mass, respec- 
tively. We assume mi = 0.1M Q and to u = 60M Q , and 
the index a is 1.35 5 . The accretion rate, equation (3), is 
therefore 

M drag ~ 1.2 x 10- 3 r; drag M bulgc (t)(l - e -w(*)), (4) 

Here we ignore the infrared luminosity from the evolved 
stars because the dust opacity for the infrared band is 
much smaller than that for the optical and UV bands. 
■Wbuige an d t(£) are directly given from the numerical sim- 
ulations. The total mass of dusty ISM accreted to the 
central massive dark object (MDO), Mmdo(*)j is obtained 

by 



Mi 



MDO 



(*) 



-J 

JO 



M diag dt. 



(5) 



As seen in equation (3), (4) and (5), the linear relation 
between the MDO mass and the bulge mass is a direct re- 
sult of the radiation drag mechanism. The possible mass 
accreted by the radiation drag in the optically thick limit 
is given by 

rti 

M M DO,max = ?7drag / L hulgc /c 2 dt ~ 5 X 10 3 Af b ulgc, 



3. RESULTS 

On the basis of the coevolution model described in the 
previous section, we estimate the mass accretion driven 
by the radiation drag during hierarchical galaxy forma- 
tion. Next, we reveal the relationship between the growth 
of a BH and that of a dark matter halo. Finally, we discuss 
the relation between AGN activity and the properties of 
the bulges. 

3.1. Mass accretion rate via radiation drag 

Evolution of the SFR and the optical depth of bulge 
(<~ 1 kpc) are shown in Figure 1. Before z <~ 4, the SFR 
and the optical depth increase, while they decrease rapidly 
for z < 4. This can be understood as follows: At high-z 
(z > 4), supply of the gas due to mergers of smaller proto- 
galaxies and consumption of the gas in the central part of 
the galaxy (bulge) are almost balanced. Thus, both the 
SFR and the optical depth are enhanced. At low-z (z < 4), 
accretion of the gas to the bulge associated with merger 
events is decreased. The gas in the bulge is consumed by 
the star formation. As a result, it makes the gas of bulge 
poor. As seen in Figure 1, we found that the evolution 
of the SFR and the optical depth are not smooth, but 
episodic. This episodic growth corresponds to the phase 
of the high mass accretion onto the bulge component trig- 
gered by major mergers with some time delays, which are 
typically 10 7 yrs (see Saitoh & Wada 2004 in details). 




(6) 

where to is 0.1 Gyr (z ~ 25) which corresponds to the 
epoch we detect the progenitor of galaxy firstly, and to is 
2.6 Gyr (z ~ 2). The correspondance between time and 
redshift is based on the cosmological model we adopted. 

In this model, we should distinguish the BH mass from 
that of an MDO, although the mass of an MDO is of- 
ten regarded as BH mass from an observational point of 
view. The radiation drag is not likely to remove the an- 
gular momentum thoroughly, and thus some residual an- 
gular momentum will terminate the radial contraction of 
the accreted gas (Sato et al. 2004). Hence, the dusty ISM 
probably forms a compact rotating torus. In this nuclear 
torus, we suppose that the mass accretion onto the BH 
horizon is determined by the Eddington rate, and that the 
BH mass grows according to 

M BH (t) = M e t / tEdd , (7) 

where iEdd is the Eddington time scale, iEdd = 
^BHMBHC 2 /iEdd, with the energy conversion efficiency, 
?7bh, and the Eddington luminosity, L^d- Unless other- 
wise stated, 77BH is assumed to be 0.42, which is the conver- 
sion efficiency of an extreme Kerr BH. Recently, Shibata 
(2004) has found that a rigidly rotating very massive star 
(VMS) with several 100M Q can be unstable for a softer 
equation of state, and eventually it forms a BH. In addi- 
tion, the theory of stellar evolution reveals that the nuclear 
burning in VMSs above 260M Q is unable to halt gravita- 
tional collapse (e.g., Heger et al. 2003). Thus, the VMSs 
inevitably evolve into massive BHs without supernova ex- ^ 
plosions. Here, we assume 260M Q as the mass of the seed ~~ r?BH 

black hole M . (8) 

5 As for the effect of IMF, If the slope and the mass range of IMF are changed to satisfy the spectrophotometric properties of galactic bulges, 
then the radiation drag efficiency is altered by a factor of ±50% (Kawakatu & Umcmura 2004). 



Fig. 1. — Redshift evolution of the star formation rate (in units 
of Mq yr _1 ) and total optical depth (rb u i ge ) of the bulge at redshift 
(z). At z > 4, both the SFR and the optical depth increase with 
time, while they decrease at low-z (z < 4). 

Figure 2 is evolution of the mass accretion rate due 
to the radiation drag (Md rag ), the rate of mass accretion 
onto a BH (Mbh)j an d the Eddington mass accretion rate 
(MEdd)- It is clear that Md rag is also episodic, reflecting 
the evolution of the SFR and optical depth (Figure 1 and 
eq.[4]). We have also found that the averaged mass accre- 
tion rate is w 10~ 5 M© yr _1 . This rate is comparable to 
the Eddington mass accretion rate for a black hole mass 
with 10 4 M n , that is , 



1 -^Edd 



Edd 



10- 5 M Q yr 



-l (mmy 1 f 

Vo.42/ I 



M BH 
10 4 M Q 
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In Figure 2, the Eddington mass accretion rate (eq. [8]) 
is larger than Md rag after z ~ 5. Since the BH mass equals 
the mass of MDO at z w 4 (see Figure 3), the mass ac- 
cretion onto the BHs after z ~ 4 would be controlled by 
the mass accretion to the MDO via the radiation drag, 
i.e. Mbh = Mdrag < Figure 2 shows that Mbh > Md rag in 
a period of 4.2 < z < 4.8 (the shaded area in Figure 2). 
In this paper, we call this period a 'BH-growing phase', 
which is w 10 8 yr. 

1.6X10 I 1 "■" ■ ~- ■ 1 1 1 1 — I 




2 4 6 8 10 12 14 

Z 

Fig. 2. — Same as Fig. 1, but for the mass accretion rate due to the 
radiation drag (Md ra g), mass accretion rate onto a BH (Mbh), and 
the Eddington mass accretion rate (Mjjdd)- The averaged mass ac- 
cretion rate due to the radiation drag is a few 10 — 5 MQyr — , which 
is comparable to the Eddington mass accretion rate for a black hole 
mass with fa 1O 4 M0. After z ~ 5, M^dd exceeds Mdrag- F° r z < 4, 
Mbh follows Mjrag because the mass of the BH reaches that of 
MDO at z ^ 4 (see also Figure 3). The shaded area (4.2 < z < 4.8) 
corresponds to the BH growing phase (Mbh > Mdrag)- 

3.2. Coevolution of MBHs, bulges and dark matter halos 

Figure 3 shows that evolution of masses of the dark mat- 
ter (Mhaio), stellar component in the bulge (Mbuige), the 
gas in the bulge (M gas ), MDO (Mmdo), and the massive 
BH (Mbh)- -Mhaio, Mb u i ge , and M gas are directly obtained 
from the numerical simulation of galaxy formation (§2.1), 

and we obtain Mmdo an d Mbh from Md rag and Mbh (Fig- 
ure 2). 

The BH mass reaches Mmdo at z « 4. As seen in Fig- 
ure 3, during z > 4, the MBH grows due to the Eddington 
mass accretion. At z < 4, the growth rate of the MDO is 
reduced because Md rag becomes smaller than Msdd, and 
-Mdrag declines owing to the decrease in the SFR and the 
optical depth of bulges (Figure 1). As a result, we find that 
an intermediate MBH with w 3 x 10 4 M Q can be formed 
at z ~ 2 in a small spiral galaxy (with a total mass of 
10 10 M Q ). At z = 2, the growth of the simulated galaxy 
finishes, so that the BH growth also would stop after z = 2 
because of the lack of gas in the bulge. In Figure 3, it is 
also found that the significantly massive dusty torus with 
Mmdo > Mbh exists in the early phase of the BH growth 
(z w 5). The massive torus itself produced by the star- 
burst in bulges would obscure the nucleus in the edge-on 
view and make a type 2 nucleus. Therefore, it is suggested 
that not only a nuclear starburst in a torus around a su- 
permassive BH like a MDO (Ohsuga & Umemura 1999; 
Wada & Norman 2002; Watabe & Umemura 2005), but 



also a starburst in a bulge may contribute to the AGN 
obscuration in the BH growing objects. 
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Fig. 3. — Same as Fig. 1, but for masses of the dark halo (Mhaio) 
, gas (Mg as ), and bulge (M^ge). Evolution of the mass of black 
hole (BH) and the massive dark object (MDO) are also plotted. 
The mass of the seed black hole is assumed to be Mo = 260Mq (see 
eq. [7]). Mb u i gc an d M ga s, respectively, while Mmdo is the mass of 
MDO and Mbh is the mass of the massive BH. It shows that the 
MDO mass is proportional to the bulge mass. The BH mass reaches 
the MDO mass at z fa 4. 
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Fig. 4. — Evolution of the mass ratios, AfBH/Mtj Ulge and 
MaH/Mhaio, against the redshift. The mass ratios are nearly con- 
stant within a factor of 2-3, but before z > 8, growth of the BH is 
slower than that of the bulge. Note that the mass ratios increase 
more rapidly during the BH growing phase (the shaded area which 
lasts fa 10 8 yr). After z ~ 4, the mass ratios are almost constant 
with M B H/M bulgc fa 5 x 10 -5 and M BH /M h ai fa 3 x 10 -6 . 

Figure 4 shows evolution of mass ratios between the 
BH and the bulge or the halo, i.e. Mra/Mbuige and 
MBH/Mhaio- We find that before z m 8, the growth of 
the bulge and dark matter halo are faster than that of 
the BH. After z m 8, the mass ratios increase with time 
because the growth of the BH is dominated by the Ed- 
dington mass accretion rate, and the rate is increasing 
exponentially (see Figure 2). One should note, however, 
that the mass ratios change by a factor of 2-3 in the BH- 
growing phase (« 10 8 yr). This change is caused by a time- 
lag between the BH growth and the growth of the bulge 
or the halo. If this is also the case in real galaxies, the 
scatter in the observed Mbh — M DU i g c relation may be be- 
cause of this time-lag. On the other hand, the mass of the 
MDO equals that of the BH [z < 4) , and the mass ratios 
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do not significantly change with MBH/-Mbuigc ~ 5 x 10~ 5 
and Mbk/M^io w 3 x 10~ 6 . Therefore, we would pre- 
dict that the scatter in the scaling relation is larger in the 
BH-growing objects at high-z than in nearby well-evolved 
galaxies. 

In Figure 5, we plot mass of the MDO and the BH 
against the halo mass. This reveals that the MDO and 
the BH coevolve with the dark matter halo from z ~ 10 
to z ~ 2. The masses of the MDO and the BH increase 
with the development of the dark halo. The mass ratio of 
the BH to the halo (Mbh/Mdm) increases gradually from 
10~ 6 to 3 x 1CT 6 from z ~ 7 to z ~ 2 (see also Figure 4). 

From these arguments, the -MBH--^haio correlation indi- 
cates that variation of the dark matter halo potential as- 
sociated with merging processes positively links with the 
mass accretion toward the galactic center via the radiation 
drag. Our model suggests that a BH mass is mutually re- 
lated to the mass of a bulge and that of a dark matter halo 
throughout the history of the galaxy formation. Compar- 
ison with observations is discussed in 84. 
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Fig. 5.— Masses of the BH (thick line) and MDO (thin line) 
are plotted against mass of the dark halo. The arrows indicate the 
masses of the dark halo at labeled redshifts. It shows that the BH 
coevolves with the dark matter halo from z ~ 10 to z ~ 2. 



3.3. AGN activity -host relation 

In this section, we examine the relation between the 
AGN activities and the properties of bulge components. 
Evolution of the bulge luminosity (itmige) a ^ optical and 
UV-band and the AGN luminosity (Lagn) are plotted in 
Figure 6. During z > 4, the AGN luminosity increases 
with the time, because the mass accretion is determined 
by the Eddington rate. After z ~ 4, the AGN luminosity is 
limited by Mdrag ( see Figure 3). Thus, the AGN luminos- 
ity exhibits a peak around z ~ 4, when Mmdo ~ Mbh- As 
seen in this figure, the AGN luminosity is always smaller 
than the bulge luminosity; in other words, no quasar 
phase, i.e. AGN luminosity dominant phase, appears. 
However, the luminosity ratio of the AGN to the bulge 
exhibits the maximal value (£agn /-^buige ^ 0.1) at z ~ 4. 
This suggests that some small spiral galaxies at high-z 
could show the same level of typical low luminosity AGNs 
in the local universe. 

Figure 7 shows the relation between the X-ray lumi- 
nosity, Lx, of the AGN and the CO luminosity, Leo, 
of the bulge. Here, the X-ray luminosity Lx is esti- 



mated assuming Lx = ex-^AGN, where ex is the X- 
ray emitting efficiency, which is supposed to be 0.1. 
The CO luminosity is derived from the gaseous mass 
assuming a conversion factor, Aco = M gas /Lco = 
4.6M Q /Kkms- 1 pc 2 (de Breuck et al. 2003). We found 
that the X-ray luminosity is positively linked with the CO 
luminosity for a wide range of luminosities, i.e. Lx ~ 
10 38 (ergs s-^Lco/Kkms" 1 kpc 2 ) 2 . 

In Figure 7, we also find that some points obviously de- 
viate from the linear relation. These points correspond to 
the BH-growing phase (4.2 < z < 4.8; the shaded area 
in Figure 2, 3, 4 and 6). This deviation implies that the 
growth of the BH via the Eddington mass accretion is 
much faster than that of gas mass in the bulge in this 
phase (see Figure 3). In other words, except for the BH- 
growing phase, our result shows that the time scale of the 
BH growth is comparable to that of the growth of the gas 
component in the bulge. Hence, by using the deviation 
from the linear relation for Lx — ^co diagram, we may 
look for the BH growing objects at high-z that we have 
missed so far. 
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Fig. 6. — AGN and bulge luminosity as a function of redshift. 
Here, we assume that Lagn is the Eddington luminosity before 
z ~ 4. After z ~ 4, the AGN luminosity is determined by the 
mass accretion due to the radiation drag. After the AGN lumi- 
nosity exhibits a peak around z ~ 4, it fades out. At z Rj 4, the 
luminosity ratio of the AGN to the bulge exhibits the maximal value 
(-^AGN/ibulgc ~ 0.1). 

4. COMPARISON WITH OBSERVATIONS 

4.1. Massive BHs in Lyman break galaxies 

In the present model, we have revealed that a small 
spiral galaxy (with a total mass of 1O 1O M0) can have an 
IMBH with ps 10 4 - 1O 5 M at z ~ 2. This is consistent 
with the mass range of the massive BH in dwarf galax- 
ies with AGNs in the local universe (Barth et al. 2004; 
Greene & Ho 2004; Barth, Greene, & Ho 2004). Moreover, 
according to Figure 6, the AGN luminosity exhibits a peak 
around z ps 4 (BH-growing phase) . This phase has several 
characteristic properties; (1) The X-ray luminosity of the 
AGN is relatively low Lx ~ 5 x 10 7 L Q if Lx = O.ILagn- 
(2) The luminosity of the bulge at optical and UV-band is 
5 x 10 9 L Q . (3) The SFR is » O.SMoyr" 1 . (4) The stel- 
lar mass is w lO 8 Af . (5) An IMBH of w 10 4 M Q exists 
in the bulge. (6) The total optical depth of the bulge is 
Tbulgo ~ 0.4. 

As for the host galaxy, it has recently been suggested 
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that Lyman break galaxies (LBGs), which are starburst 
galaxies at high-z, may be in the forming phase of a galac- 
tic bulge (Friaca & Terlevich 1999; Matteucci & Pipino 
2002). The SFR in the LBGs is - 3 - SOOMoyr" 1 . LBGs 
have a typical luminosity ~ 10 10 — 10 11 L Q , a stellar mass 
~ 10 10 — lO n M , and they are observed to be optically 
thin (e.g., Shapley et al. 2001). They also exhibit strong 
clustering at z^ 3, suggesting that hierarchical clustering 
is on-going. In addition, the Chandra X-ray observatory 
has detected the hard X-ray of LBGs at z = 2 — 4 with 
the luminosity of <~ 1O 8 L although it is still uncertain 
whether the X-ray emission arises from AGNs (Brandt et 
al. 2001). Comparing these properties of LBGs with our 
predictions (i.e. (l)-(6)), the small spiral galaxy (with a 
total mass of 1O 1O M0) in the BH-growing phase may cor- 
respond to low-mass counterparts of the LBGs. Extrapo- 
lating the scaling relation that we found to the observed 
LBGs, they would have the MBHs with w 10 6 - 10 7 M Q . 
According to recent observations, only 3% of LBGs show 
AGN activity in the rest-frame hard X-ray band (Nam- 
dra et al. 2002) and optical band (Shapley et al. 2001). 
However, it has not been clear that this 3% fraction re- 
flects the duty cycle of mass accretion to BHs (see Shap- 
ley et al. 2001 in details) or the possibility that the LBGs 
have massive BHs. Hosokawa (2004) claimed, assuming 
MeH/Mbuigc ~ 0.001, that ~ 10% of LBGs at z ~ 3 can 
have a massive BH with <~ 1O 7 M to reproduce the lo- 
cal mass function of SMBHs (Salucci et al. 1999; Yu & 
Tremaine 2002; Aller & Richstone 2002; Shankar et al. 
2004). 

4.2. Mbh — A^buige ond Mbh — Mhaio relations 

Barth et al. (2004) suggested that the correlation 
between BH mass and stellar velocity dispersion (the 
Mbh — o~ relation) holds on a mass scale of an intermediate 
BH with w 10 4 - 10 6 M Q . In addition, Baes et al. (2003) 
found a correlation between the BH mass and halo mass, 
i.e. M BH /1O 8 M ~ O.lltMhaio/lO^Mo) 1 - 27 , by using the 
Mbh - cr relation. This relation gives M B h = 3 x 1O 4 M 
for Mhaio = 1O 1O M , which is comparable to our predic- 
tion (Mbh ~ 3 x 1O 4 M at z <~ 2). However, it should be 
noted that there are large uncertainties in the empirical 
laws, as mentioned by Ferrarese (2002). 

In our model, the final BH mass-to-bulgc mass ratio 
is w 5 x 10~ 5 , which is much smaller than the observed 
value (ps 10~ 3 ) in nearby large galaxies. Thus, if our sce- 
natio is correct, it is expected that small spiral galaxies 
at high-z, which are not direct counter parts of the dwarf 
galaxies at low-z, could have the IMBHs and the smaller 
BH mass-to-bulge mass ratios. Currently, it is difficult to 
detect the IMBHs in small galaxies at high-z, and there- 
fore we can not directly prove our prediction, namely the 
small BH mass-to-bulge mass ratio. Moreover, there is 
also a room in the theoretical model, especially on effects 
of mechanical, radiative and chemical feedback processes 
from star formation. For instance, it is not trivial whether 
the stellar feedback affects on the SFR positively or neg- 
atively. Therefore, it is ultimately necessary to perform 
high-resolution radiative hydrodynamical simulations for 
galaxy formation taking into account these effects explic- 
itly. " 

In this paper, we focus on the formation of a MBH due 



to the radiation drag in a small spiral galaxy. However, the 
process discussed here could be applied to more massive 
galaxies because the transfer of the angular momentum 
via the radiation drag is independent of the mass scale of 
the galaxies. By using the BH mass-to-halo mass relation 
for a small spiral galaxy (MsH/Mhaio w 10~ 6 ), we can pre- 
dict that the massive spiral galaxies with Mh a i = 1O 12 M 
have MBHs with 10 6 - 1O 7 M . In addition, the observa- 
tions have suggested the SFRs of the massive galaxies at 
high-z are 10-100 times higher than those of small galaxies 
like the one considered here. Thus, the final BHs might 
achieve w 1O 8 M because the effect of the radiation drag 
is linearly proportional to the SFR (eq.[4]). If this is the 
case, host galaxies of the luminous quasars at high red- 
shift would be star-forming or post-starburst galaxies. On 
the other hand, we should note that the mass of a BH 
would depend on the morphology of the galaxies even if 
the mass of the dark halo is the same. This is because 
the rate of mass accretion onto the BHs is not determined 
by the disk components, but by the bulge components in 
the host galaxies, due to the effects of geometrical dilu- 
tion and opacity (see Kawakatu & Umemura 2004 in de- 
tails). Therefore, the morphology differences of host galax- 
ies would cause the large scatter in the mass relation of the 
BHs to the halos in spiral galaxies. In fact, some authors 
claim that the Mbh — Mhaio relation is much weaker than 
the Mbh — Mb u i ge relation in spiral galaxies (Salucci et al. 
2000; Zasov et al. 2004). 
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Fig. 7. — The relation between the CO luminosity (Leo) and 
the X-ray luminosity (Lx) during the evolution of the BH and 
the galaxy. The filled circles show time evolution in our model. 
The X-ray luminosity is positively correlated with the CO luminos- 
ity. The thin line represents this correlation, i.e. logLx(crg/s) ~ 
38 + 2 log Leo (K km s ~ 1 Kpc 2 ) . The shaded area corresponds to the 
BH growing phase (4.2 < z < 4.8). The arrows show time evolution 
in the BH growing phase. 

4.3. Lx — Leo relation 

In §3.3, we predicted that the X-ray luminosity of AGNs 
is positively correlated with the CO luminosity of bulges 
from z ~ 10 to z ~ 2. This correlation (i.e. Lx ~ 
10 38 (ergs s" 1 ) 

(Lco/Kkms -1 kpc 2 ) 2 ), if we extrapolate it to AGNs with 
higher luminosity, is consistent with the Lx — Leo rela- 
tion found in the low rcdshift Seyfert galaxies and quasars 
(Yamada 1994). The observed Lx — Lqo relation shows 
a scatter of about one order of magnitude. Suppose all 
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galaxies follow the same Lx — -^co relation that we found 
here, we suspect that the large scatter in the observed 
scaling relation was caused when the BHs were in their 
growing phases. Submillimeter observations by ALMA for 
luminous high-z quasars will be helpful to investigate the 
connection between the black hole mass and the host prop- 
erties. 

5. CONCLUSIONS 

Combining a theoretical model of the mass accretion 
onto a galactic center due to the radiation drag with high- 
resolution 7V-body/SPH simulations (2 x 10 6 particles, one 
SPH particle has 10 3 M Q and a softening length of <~ 50 
pc), we demonstrate growth and formation of a massive 
BH during hierarchical formation of a small spiral galaxy 
(with a total mass of 10 10 M Q ). We found that the aver- 
age rate of the mass accretion due to the radiation drag 
is w 10 _5 M Q yr _1 . Finally, a small spiral galaxy can have 
an IMBH with 10 4 - 1O 5 M at z ~ 4. 

Our model suggests that the growth of the massive BHs 
correlates not only with that of the galactic bulges, but 
also with that of the dark matter halos in the hierarchical 
formation of spiral galaxies. The massive BHs coevolve 
with the dark matter halo from z^l5toz^2. This 
means that the change in the dark matter potential closely 
correlates with the rate of the mass accretion onto a seed 
BH with the help of the radiation drag. The final mass ra- 
tio of the BH-to-dark matter halo is s=s 10~ 6 and the final 



BH-to-bulge mass ratio is about 5 x 10~ 5 in a small spi- 
ral galaxy, which is much smaller than the observed value 
(w 10~ 3 ) in the large galaxies due to the opacity effect, 
although the stellar feedback would affect on the result. 
Moreover, the time-lag between the BH growth and the 
growth of the bulge (halo) would cause the scatter of the 
observed scaling relation. 

In terms of relationship between the AGN activity and 
the properties of host galaxies, we found that even a small 
spiral galaxy could show the same level of the typical low 
luminosity AGNs with LAGN/^buigc w 0.1 at z 4. Our 
model shows that the X-ray luminosity of the AGN is 
positively correlated with the CO luminosity (the gaseous 
mass) of the bulge very well. Furthermore, our result pre- 
dicts that the BH-growing objects deviate from this scaling 
relation. By comparing our results with the properties of 
the LBGs, we predict that the LBGs could harbor massive 
BHs with 10 6 - 10 7 M n . 
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